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I.  INTRODUCTION 


SCOPE  OF  THE  CONTRACT 

1.1  The  purpose  of  this  contract  is  the  study  and  development  of 
an  angular  rate  indicating  instrument  based  on  measuring  the  angular 
rate  of  an  artificially  induced  vortex.  The  measurement  of  the  change 
in  rate  of  the  artificially  induced  vortex  is  the  measure  of  an  externally 
superimposed  rotation  on  the  system.  The  successful  achievement 

of  this  development  would  provide  an  instrument  to  replace  gyroscopes. 

PHASES  OF  THE  CONTRACT 

1.2  The  scope  of  the  work  in  Contract  NOw  62  -  0770  -  c  was 
written  in  three  phases  to  insure  the  orderly  conduct  of  the  project. 
They  read  as  follows: 

Phase  A  -  Conduct  investigations  limited  to  veri*  'ation  of 
practical  methods  of  measuring  the  rate  of  change  of  gas  ve  cities 
necessary  to  successful  operation  of  a  vortex  gyroscope. 

Phase  B  -  Conduct  investigations  limited  to  verification  of 
practical  methods  of  providing  constant  initial  gas  parameters  neces¬ 
sary  to  successful  operation  of  a  vortex  gyroscope. 

Phase  C  -  Complete  work  proposed  in  Development 
Laboratories  proposals  Nl.  The  contractor  shall  limit  research  and 
development  effort  initially  to  Phase  A  until  given  written  instructions 
by  the  technical  director  to  proceed  to  Phase  B.  Similar  instructions 
shall  be  obtained  before  proceeding  from  Phase  B.  to  Phase  C. 
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II.  WORK  ACCOMPLISHED  DURING  FIRST  QUARTER 


ORGANIZATION  OF  THE  WORK 

2.1  Activities  on  Contract  NOw  62  -  0770  -  c  started  the  15th 
of  March,  1962.  This  report  covers  the  work  accomplished  by  the 
15th  of  June,  19  62.  Figure  1  shows  the  organization  of  the  work 
in  bar  chart  form.  The  work  carried  out  during  the  first  quarter 
followed  this  chart  closely. 

STUDIES  ON  MIXED  ION  SYSTEMS 

2.2  The  studies  carried  out  on  mixed  ion  systems  are  presented 
in  Appendix  A.  The  studies  performed  to  date  have  lead  to  the 
following  conclusions: 

The  heavy  ion  diffusion  rates  dictate  an  environ¬ 
ment  near  car  above  normal  gas  density.  At  signifi¬ 
cantly  lower  density,  the  diffusion  rates  become 
large  and  the  pulse  or  wave  shape  of  the  group  is 
rapidly  distorted. 

If  mixed  ion  systems  (heavy  positive  ions  plus 
electrons)  are  considered,  then  recombination  rates 
at  normal  gas  density  are  high  and  it  is  difficult, 
maybe  impossible,  to  retain  a  measurable  ion  group 
at  the  down-stream  observation  point.  This  obser¬ 
vation,  derived  from  studies  of  Xe  ions  in  a  gas  of  normal 
density,  is  less  true  of  a  helium  system  at  10  atmospheres. 

It  cannot  be  concluded  from  the  work  done  that  systems 
with  mixed  ions  cannot  be  made  to  work.  However, 
more  promising  concepts  have  arisen  which  are  not 
subject  to  the  recombination  limitation.  The  exploration 
of  these  concepts  appears  to  be  more  fruitful  than  a 
continued  study  of  the  mixed  ion  systems. 
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STUDY  ON  SINGLE  ION  SYSTEMS 


2.3  The  studies  carried  out  on  the  diffusion  of  single  ion  systems  are 
presented  in  Appendix  B.  The  studies  performed  to  date  have  lead  to  the 
general  conclusion  that  the  single  ion  system  does  not  have  the  inherent 
limitation  imposed  by  recombination  in  the  mixed  ion  case,  and  that  a 
system  based  on  this  could  be  successfully  developed. 

STUDIES  ON  THERMAL  SYSTEMS 

2.4  The  studies  carried  out  on  the  thermal  system  are  presented  in 
Appendix  C.  The  thermal  system  consists  of  pulsed  heater  and  a  down¬ 
stream  detector  which  senses  the  arrival  of  pulses  of  heated  gas.  The 
detector  may  be  a  thermocouple  or  a  resistance  thermometer.  The  heater 
can  be  a  simple  thin  wire.  The  studies  performed  have  lead  to  the 
general  conclusion  that  this  system  should  be  far  the  easiest  of  devel¬ 
opment  and  that  it  can  be  made  to  work  adequately. 

THE  GAS  VELOCITY  MEASURING  SETUP 

2.5  A  test  setup  was  designed  and  constructed  during  the  three  months 
period  being  reported.  This  setup  consists  of  equipment  developed  for 
regulating  the  passage  of  the  working  fluid  in  a  test  section  wherein  various 
electrode  systems  can  be  placed  for  testing.  This  test  setup  is  being  used 
in  conjunction  with  a  cathode  ray  oscilloscope  and  associated  electronic 
components  for  the  visualization  of  the  signals  transmitted  and  received  by 
the  electrode  system  under  test. 

TEST  RESULTS 

2.6  After  the  usual  initial  difficulties  in  getting  a  system  functioning, 
useful  results  on  the  measurement  on  the  working  fluid  velocity  are  being 
obtained.  These  results  are  being  used  to  confirm  theoretical  expected 
values  and  establish  threshold  of  operation  of  measuring  electrode  systems 
with  various  working  fluids.  These  tests  will  continue  during  the  next 
quarter . 
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III.  THE  WORK  PLAN 


THE  WORK  PLAN 

3.1  Since  the  problem  associated  with  Phase  A  of  the  work 
statement  is  being  satisfactorily  solved  and  since  the  working 
fluid  pump  problem  has  been  circumvented  by  a  mechanical 
pump  in  place  of  the  original  method  suggested  in  our  proposal, 
the  work  plan  has  been  modified  to  the  one  shown  in  Figure  2, 
which  is  self-explanatory. 

EXPECTED  WORK  TO  BE  ACCOMPLISHED  DURING  NEXT  QUARTER 

3.2  By  the  end  of  the  next  quarter,  we  hope  to  accomplish  most 
of  the  theoretical  work  associated  with  the  signal  generation  and 
detection  and  be  well  on  our  way  into  the  motion  analysis  study 
of  a  complete  angular  rate  indicator,  as  well  as  having  considered 
various  alternatives  relating  to  the  electrode  geometry  and  space. 
We  will  hope  to  have  finished  all  flow  velocity  measurements  that 
are  required  and  that  the  testing  of  an  angular  rate  instrument  be 
under  way.  This  of  course  means  that  the  design  of  an  instrument 
be  also  completed  by  then. 
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Design  of  basic  A.R.I.  instrument  for  test 
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Figure 


APPENDIX  A 


STUDIES  PERFORMED  ON  MIXED  ION  SYSTEMS 

1.1  The  Development  Laboratories  Proposal  N-l  suggested  that  gas 
speeds  might  be  accurately  measured  by  injecting  or  creating  low  speed 
ions  into  the  gas  stream  and  observing  the  arrival  time  of  these  ions  at 
some  down-stream  point.  The  studies  performed  to  date  have  led  to  the 
following  conclusions: 

1.  The  heavy  ion  diffusion  rates  dictate  an 
environment  near  or  above  normal  gas  density.  At 
significantly  lower  density,  the  diffusion  rates 
become  large  and  the  pulse  or  wave  shape  of  the 
group  is  rapidly  distorted. 

2.  If  mixed  ion  systems  (heavy  positive  ions  plus 
electrons)  are  considered,  then  recombination  rates  at 
normal  gas  density  are  high  and  it  is  difficult,  maybe 
impossible,  to  retain  a  measurable  ion  group  at  the 
down-stream  observation  point.  This  observation, 
derived  from  studies  of  Xe  ions  in  a  gas  of  normal 
density,  is  less  true  of  a  helium  system  at  10  atmospheres. 

3.  Detection  efficiencies  have  not  been  explored 
and  raise  a  major  uncertainty  in  using  the  above  con¬ 
clusions  to  make  decisions.  Detection  systems  will  be 
somewhat  different  in  mixed  ion  versus  single  ion  systems. 

Future  analytical  work  should  endeavor  to  establish  the 
detection  efficiency  in  single  ion  systems. 

4.  The  concurrent  development  of  a  system  using 
heated  gas  pulses  detected  with  a  thermocouple  is  already 
in  the  experimental  stage.  At  the  present  time,  this  work 
is  showing  sufficient  promise  to  favor  its  support  instead 
of  the  originally  proposed  ion  scheme. 

1.2  The  remainder  of  Appendix  A  and  also  Appendix  B  describe  the  theoretical 
analysis  which  have  been  performed  upon  ion  systems. 
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MINIMUM  MEASURABLE  CURRENTS 


1.3  The  geometry  and  gas  rates  are  those  described  in  Development 
Laboratories  proposal  N-l.  A  gas  speed  of  about  15  cm/second  is  to  be 
observed  over  a  transit  path  of  about  0.5  cm  by  injecting  ions  and 
measuring  their  speed  as  they  drift  with  the  gas.  The  pulses  then 
arrive  at  the  observation  point  at  33  millisecond  intervals  yielding  a 
fundamental  frequency  of  30  cps.  If  we  assume  an  amplifier  input  pass 
band  of  10  to  100  cps  and  an  input  resistance  of  5  megohms,  then  the 
mean  square  noise  current  is:^ 

-TZ  =  4  k  T  .  B  (1) 

R 

B  =90  cps 
6 

R  =  5  x  10  chms 

kT  =  4.14  x  10  joules  at  300°  k 

—13 

Accordingly,  i  (rms)  =  5.5  x  10  amperes.  The  corresponding  voltage 
across  the  5  megohm  input  resistance  is  2.75  microvolts.  It  may  be 
noted  that  although  reduction  of  the  input  band  pass  to  0.9  cps  will 
reduce  the  noise  current  ten-fold,  that  one  is  thereby  inviting  electronic 
refinements  which  are  not  entirely  within  the  spirit  of  the  simple  gyro 
substitute . 

1.4  The  ion  detector  is  regarded  as  a  pair  of  parallel  electrodes  1  cm 
long  of  very  thin  wire  and  spaced  at  0.05  cm  apart.  It  will  be  assumed 
that  the  collection  aperture  is  1  cm  by  0.05  cm  =  0.05  cm  .  The 
detector  is  instrumented  to  sense  changes  in  the  conductivity  of  the 
gas  and  it  is  assumed  that  all  of  the  positive  ions  in  the  aperture  are 
collected  at  one  electrode  and  all  of  the  electrons  are  collected  at  the 
other  electrode.  The  total  current  through  the  electrode  systems  is  then: 

i/  J.J.  Freeman,  Principles  of  Noise,  Wiley  and  Sons,  1958, 

p.  115. 
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1.6  In  progressing  from  the  point  of  creation  to  the  detector,  ion 
densities  are  reduced  by  diffusion  and  by  recombination.  Whereas 
diffusional  losses  are  rather  easily  controlled,  recombination  rates 
pose  a  serious  problem  to  sustaining  the  minimum  detectable  current 
in  the  above  geometry. 

ION  DIFFUSION  (WITHOUT  RECOMBINATION) 

1.7  It  can  be  shown  that  the  electrons  cannot  diffuse  away  from 
the  ions  because  the  electrical  potential  energies  would  then  greatly 
exceed  the  available  kinetic  thermal  energy.  There  is  some  distortion 
to  the  diffusion  solution  in  this  system,  however  we  shall  first  treat 
the  ions  as  if  electrical  forces  were  negligible.  A  simple  model  is 
conceived  in  which  ions  and  electrons  are  created  in  equal  number  in 
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a  nearly  line  source  geometry  which  then  is  dragged  into  a  sheet  by  the 
passing  gas.  This  is  illustrated  in  Figure  1  in  which  the  gas  velocity 
is  parallel  to  the  time  axis,  the  vertical  axis  represents  -Ke  ion  density. 


FIGURE  I. 

The  initial  width  of  the  sheet  is  taken  as  0.01  cm  (2%  of  the  transit  path 
length).  The  ion  density  is  modulated  sinusoidally  at  about  30  cps.  The 
ion  source  is  conceived  as  a  narrow  pencil  (0.01  cm  diameter)  of  UV  light 
modulated  by  a  graded  disc  driven  at  1800  rpm. 

1.8  As  the  sheet  drifts  from  its  point  of  creation  to  the  detection  point,  it 

spreads  laterally  by  diffusion.  The  diffusion  constant,  at  normal  density 

(1  atmosphere) ,  for  positive  xenon  ions  in  xenon  is  estimated  from  mobility 

2/ 

data  to  be: 

D  =  1.35  xlO  2  cm2  /sec 


The  mobility  data  is  that  of  R.N.  Varney  as  reported  by  S.C.  Brown, 
Basic  Data  of  Plasma  Physics.  Wiley  and  Sons ,  1959,  0.  72, 

Figure  3.35.  The  value  of  D  given  abov§^s  consistent  with  D=  \v/3 
and  a  collision  cross-section  of  2.1  x  10  cm  yatom. 
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If  one  assumes  a  simple  solution  of  the  plane  source  diffusion  equation 
then  the  ions  are  contained  between  parallel  planes  at  x  =  +  V  4Dt. 
Assuming  one  square  centimeter  of  area  (normal  to  the  diffusion  direction, 
x) ,  the  system  volume  at  time  t  is: 

V  =  2  x  =  2  V  4  Dt  (3) 


If  this  volume  contains  N  particles  uniformly  distributed  and  if  no 
particles  disappear  by  recombination,  the  particle  density  at  time  t  is: 

N  _  N 

n  - - 

2  V  4  Dt  (4) 


The  singularity  at  t  =  0  is  characteristic  of  singular  source  diffusion 

problems  and  is  removed  by  computing  the  time  t  for  the  plane  source 

to  diffuse  to  the  initial  finite  thickness  (2  x  =0.01  cm).  From  (3)» 

o 

this  time  is  found  to  be  0.463  milliseconds.  Accordingly,  Equation  (4) 
will  be  rewritten: 


where  n  =  the  initial  concentration  at  creation, 
o 


1.9  Table  I  shows  densities  33.3  milliseconds  after  creation  (0.5  cm 
downstream  from  the  creation  point)  as  a  function  of  initial  density  and 
operating  pressure.  Operating  pressure  influences  the  diffusion  constant 


and  hence  t  in  accordance  with 
o 

t'  -  -fi-ia 

o  D' 


=  -Pjto.  ■ 
p 
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TABLE  I 


DENSITY  REDUCTION  BY  DIFFUSION  ONLY 


Ion  Pairs/cm°  after  33.3  milliseconds 
p  =  1/2  atm _ p  =  1  atm _ p  =  2  atm 


n  (ions/cm8 
o _ 

5  x  107 
5  x  108 

9 

5  x  10 


3.45  x  105 

3.45  x  106 

3.45  x  107 


6.85  x  105 

6.85  x  106 
7 

6.85  x  10 


1.35  x 106 

1.35  x 107 

1.35  x 108 


It  may  be  noted  that  (a)  the  density  reduction  is  strongly  Influenced  and 
easily  controlled  by  choice  of  an  operating  pressure;  (b)  that  final 
density  is  directly  proportional  to  initial  density;  and  (c)  that  the  mini- 

g 

mum  observable  density  of  1.53  x  10  ions  is  met  with  most  of  the  initial 
densities  and  pressures  of  Table  1. 

1.10  When  recombinations  are  taken  into  account,  it  is  found  that 
final  density  is  no  longer  proportional  to  initial  density  and  that  it 
becomes  more  difficult  to  sustain  the  minimum  observable  density. 

DIFFUSION  WITH  RECOMBINATION 

1.11  If  there  are  N  total  particles  in  volume,  V  then  the  density  n  =  N/V. 
The  density  change  with  time  is  given  by: 

_dn_  _  _  _£L_  aV  ,  1  ( ±  N)  (7) 

dt  v2"  ^ t  v  Vat/ 

N  /v 

The  partial  derivatives  are  taken  at  constant  N  and  constant  V  respectively. 
With  equation  (3)  and  n  =  N/V,  this  equation  can  be  written: 


dn  _  _n_  1  l  d  N  \  , 

dt  ”  '  t  + 1  +  v  \at  ) 

with  oN/  d  t  =0,  this  equation  has  the  solution  (5).  The  partial 
^  N/  ^  t  at  constant  volume  is  the  particle  loss  due  to  recombination  . 
It  may  be  expressed  by: 


AM 

dt 


=  -  oC-V 


The  recombination  constant, 


is  taken  as 


2  x  10  6cm3/sec.^// 


(9) 


Using  this  equation  (9)  becomes: 


dn 

dt 


n 

t  +  t 

o 


2 

-  ccn 


(10) 


As  before,  the  symbol  tQ  represents  the  time  required  for  a  plane  source 
to  diffuse  to  its  initial,  finite  half  thickness,  0.005  cm,  and  is  0.463 
milliseconds  for  operation  at  o  atmosphere. 


1.12  The  solution  of  (4)  if  n  is  the  density  at  which  particles  are  created  is: 

o 


1 

n 


+  oc  t  In 
o 


(ID 


The  following  table  shows  the  particle  density  at  t  = 
milliseconds  as  a  function  of  n  when  oc  =  2  x  10 


milliseconds . 


o 


10,  20, 

and  t 

o 


and  33.3 
=  0.463 


S.C.  Brown,  loc.  cit,  p. 


195,  Table  8.3. 


see  p.  188 


For  the  source  of  equation  (9) 
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TABLE  II 


RESIDUAL  ION  DENSITIES  AT  DOWN  STREAM  OBSERVATION  POINTS 


n 

o 

3 

Ion  Pairs  Per  Cm 

t  =  10  ms 

t  =  20  ms 

t  =  33 

.3  ms 

5  x  107 

1.93  x  106 

0.98  x  106 

0.57 

x  106 

5  x  108 

9.05  x  106 

4.10  xlO6 

2.30 

x  106 

9 

5  x  10 

1.43  x 107 

6. 24  x  106 

3.29 

x  106 

oo 

1.53  x 107 

6.45  x  106 

3.45 

x  106 

1.13  It  was  previously  shown  that  the  particle  density  at  the  observation 
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point  should  be  1.53  x  10  ion  pairs/cm  .  The  table  above  shows  that 
although  this  can  be  attained,  there  is  not  much  leeway  if  times  as  long 
as  33.3  milliseconds  are  involved.  Increasing  the  pressure  appears  to 
be  slightly  advantageous ,  this  is  shown  in  Table  3  which  applies  only  to 
observations  made  at  33.3  milliseconds.  At  the  higher  pressures,  recombi¬ 
nation  is  the  dominant  influence ,  and  at  nQ  =  00  ,  recombination  is  the 
only  factor  responsible  for  reducing  particle  density  with  time. 


TABLE  III 


EFFECT  OF  PRESSURE  UPON  RESIDUAL  ION  PAIR  DENSITY  AT  33.3  MILLISECONDS 


3 

Ion  Pairs/cm 

after  33.3  Milliseconds 

p  (atm) 

2 

D  (cm  /sec 

n  -  5  x  10  2 

o 

n  =  5  x  10^ 
o 

n  = 
o 

9 

5  x  10 

1/2 

2.70  xlO-2 

3.10  x  105 

1.60  x  106 

2.75 

x  106 

l 

1.35  x  10“2 

5.72  x  105 

2.30  x  106 

3.29 

x  106 

2 

6.75  x  10"3 

1.01  x  106 

3.11  x  106 

3.93 

x  106 

5 

2.70  x  10"3 

1.98  x  106 

4.44  x  106 

5.05 

x  106 

10 

1.35  x  10'3 

3.08  xlO6 

5.68  x  106 

6.19 

x  106 

50 

2.70  x  10~4 

6.91  x  106 

9.47  x  106 

9.88 

x  106 

oo 

0 

1.15  x  107 

1.46  x  107 

1.50 

x  107 

1.14  Of  particular  note  is  the  fact  that  at  initial  densities  greater  than 
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5  x  10  per  cm  ,  the  residual  density  at  33.3  milliseconds  is  very  unrespon¬ 
sive  to  further  increases  in  initial  density  (compare  the  last  two  columns 
of  Table  3).  For  instance,  a  computation  with  nQ  =  5  x  10^  yields  essentially 

the  same  result  as  that  at  n  =  5  x  10^  . 

o 

CONCLUSIONS 

Because  recombinations  play  a  dominant  role  in  limiting  the  final 
ion  density,  the  choice  of  xenon  was  perhaps  unfortunate.  Helium  with 
a  recombination  rate  two  orders  of  magnitude  lower  than  xenon,  operated 

at  10  acmospheres  to  limit  diffusion  rates,  would  reduce  in  density  from 

9  7  4/ 

5  x  10  to  about  9 . 5  x  10  at  33 . 3  ms . 


4/ 


Assuming  the  diffusion  constant  is  inversely  proportional  to  the  square 
root  of  the  mass  and  inversely  proportional  to  the  operating  pressure. 
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It  cannot  be  concluded  from  the  foregoing  that  systems  with 
mixed  ions  cannot  be  made  to  work.  However,  more  promising  concepts 
have  arisen  which  are  not  subject  to  the  recombination  limitation. 

The  exploration  of  these  concepts  appears  to  be  more  fruitful  than  a 
continued  study  of  the  mixed  ion  systems. 
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APPENDIX  B 


STUDIES  PERFORMED  ON  SINGLE  ION  SYSTEMS 

1.  The  studies  made  of  mixed  ion  systems  indicated  that  ion 
recombinations  would  occur  rapidly  and  that  the  residual  ion  density 
at  the  point  of  observation  might  be  only  marginally  observable.  If 

ions  of  a  single  sign  are  considered  (say  electrons  captured  by  heavy  ions 
such  as  bromine)  then  recombination  problems  are  disposed  of.  One 
implementation  of  this  scheme  involves  emitting  electrons  from  a 
thoriated  tungsten  filament  into  a  gas  having  a  high  electron  attachment 
probability.  Although  not  extensively  examined ,  this  method  is  regarded 
as  a  more  fruitful  line  of  endeavor  than  the  mixed  ion  systems. 

2.  In  the  previously  considered  mixed  ion  systems,  the  conventional 
picture  of  diffusion  was  preserved  because  of  the  tendency  to  microscopic 
electrical  neutrality  throughout  the  system.  This  is  not  the  case  with 
ions  of  a  single  sign  and ,  at  least  under  some  conditions ,  the  effect  of 
electrical  forces  can  greatly  exceed  the  diffusion  effects.  The  analysis 
of  this  problem  proved  to  be  rather  lengthy  and  is  briefly  outlined  at  the 
end  of  this  Appendix  (Note  1). 

3.  A  system  is  conceived  in  which  electrons  are  emitted  in  a  short 

burst  from  a  heated  filament  and  are  almost  instantly  captured  by  a  gas 

having  high  electron  affinity.  The  negative  ions  so  formed  then  drift 

with  the  main  gas  stream  to  a  detection  grid  where  their  arrival  is  sensed 

by  the  passage  of  a  current  in  an  external  circuit.  The  minimum  detectable 

current  is  similar  to  that  computed  for  the  mixed  ion  system,  i.e.  5.5  x  10  * 

-19 

amperes.  With  the  electron  charge  equal  to  4.8  x  10  coulombs,  the 
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number  of  ions  emitted  by  the  source  ist 


N 


5.5  x  10  ^  !_  = 

4.8  x  10~19  X  30 


3.  82  x  10 


4 


per  cycle  per  cm  of 
filament  length. 


The  ions  are  formed  rapidly  and  diffuse  cylindrically  from  the  region 
around  the  filament  (See  Note  2).  The  pulse  rate  is  assumed  to  be  30 
cycles  per  second. 


4. 


It  is  shown  in  Note  1  that  the  quantity  N6  /kT  may  be  used 


as  a  measure  of  the  relative  importance  of  space  charge  effects  and 

diffusion  effects  in  the  expansion  of  the  ion  cylinder.  At  300°  K  , 
-14  2  -20 

k  T  =  4.14  x  10  ergs,  €  =  2.56  x  10  erg-cm.  Hence 


N  6  /  k  T  = 


-4  -70 

(3.82  x  10  )(2.56  x  10  ) 

4.14  x  10-14 


=  0.0236 


Using  Equation  9  of  Note  1,  it  is  apparent  that  the  electrical  forces 
are  slight  at  this  initial  concentration  and  that  the  system  can  be 
described  by  the  usual  diffusion  relations.  The  initial  concentration 
could  be  increased  10  fold  without  seriously  distorting  the  diffusion 
viewpoint . 


5  No  data  has  been  found  describing  the  diffusion  and/or 

mobility  of  halogen  ions.  Data  on  0  in  09  might  be  expected  to  be 

2  9/ 

somewhat  similar  to  fluorine  and  are  available.  For  0  in  0  ,  the 
mobility  is  about  2.5  cm/sec  per  volt/cm  at  1  atmosphere  and  273°  K. 
The  diffusion  constant  and  mobility  are  simply  related  by: 


A 


.,e_ LL 
k  T 


(10) 


2/ 

3/ 


S.C.  Brown,  "Basic  Data  of  Plasma  Physics , "  Wiley  and  Sons,  1959, 
pp.  68-69 

? 

Consistent  units  for  this  equation  in  the  cgs  system  /-o,  cm  /sec- 
statvolt;  e,  1.6  x  10  10  esu;  k,  ergs  per  degree  Kelvin;  T  in  °K;  and 
D  in  cm2/sec. 
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This  yields  D  =  0.176  cm  /sec  for  0  in  0  and  is  taken  as  an  estimate 

Zi 

of  the  diffusion  constant  of  F  in  F^.  Fluorine  is  of  interest  because  of  its 
high  election  affinity.  It  is  probable  that  a  small  amount  of  a  heavier 
halogen  in  a  heavy  inert  gas  such  as  xenon  would  provide  a  lower,  more 
desirable,  diffusion  constant.  In  view  of  these  comments ,  it  is  reasonable 
to  believe  that  a  diffusion  constant  of  0.05  cm  /sec  can  be  achieved  without 
increasing  the  system  pressure  above  1  atmosphere.  With  such  a  diffusion 
constant,  the  cylinder  would  expand  from  a  line  source  to  about  0.08/cm 
radius  during  its  33  millisecond  transit  time  to  the  detector  and  its  arrival 
at  the  detection  grid  would  occur  over  a  time  interval  of  about  10  milliseconds. 
The  detection  grid  should  be  on  the  order  of  0.16  cm  wide  in  order  to  capture 
a  substantial  portion  of  the  negative  ions  transmitted  from  the  source. 

6.  Because  of  the  finite  distance  the  electrons  travel  prior  to 
attachment,  the  initial  configuration  of  diffusing,  negative  ions,  is  finite. 
However,  the  rough  calculation  of  Note  2  indicates  that  a  few  per  cent  of 
halogen  in  xenon  will  attach  90%  of  the  initial  electrons  in  a  radius  which 
is  smaller  than  the  diffusion  distance  (0.08  cm)  estimated  above. 

7 .  The  foregoing  scheme  using  negative  ions  appears  to  be  flexible, 
composition  and  pressure  can  be  varied  to  control  diffusion,  and  the 
scheme  does  not  have  the  inherent  limitation  imposed  by  recombination  in 
the  mixed  ion  case. 
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NOTE  1: 


The  Effect  of  Space  Charge  Upon  Diffusion  Rate 


A  general  equation  for  diffusion  with  forces  was  derived  by 

M.H.  Johnson^.  With  the  aid  of  hydrodynamic  continuity,  assuming 

isothermal  (p  =  nkT)  behavior  and  irrotational  flow,  and  neglecting  the 

2  2 

streaming  kinetic  energy  and  b  n/  d  t  ,  it  is  possible  to  reduce  this 
equation  in  the  case  of  cylindrical  symmetry  to 


1 

r 


A 

3r 


(r 


,  r  n  F  '  _J_  ^n 

dr  V  kT  }  D  at 


(1) 


Here  n  is  the  particle  density  at  radius  r  and  time  t,  F  is  the  force  acting 
on  each  particle,  D  is  the  diffusion  constant,  T  the  temperature  and  k  is 
Boltzmann's  constant. 


The  force  acting  on  a  particle  at  r  is  due  to  the  charge  interior 
to  r  and  may  be  written: 


F  = 


4  ire  2 
r 


f  r 
J  o 


n(r')  r'  dr' 


(2  ) 


where  6  is  the  electronic  charge  (it  is  assumed  that  each  charged  particle 
is  only  singly  charged). 


Inserting  Equation  (2)  in  (1)  the  equation  for  diffusion  with  space 
charge  is: 

1  ^  3  n  >  _  _1_  _jd_ 

r  ^r  fcir  r  £  r 


ATT  G 
k  T 


»f  n 
J 


(r1)  r*  dr' 


+  -Jr  — n  (3) 
D  at 


See  the  summary  in  Chapman  and  Cowling,  Mathematical  Theory  of 
Non-Uniform  Gases,  Second  Edition,  Cambridge  University  Press, 
p.  415,  Equation  1. 
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It  will  be  seen  that  if  the  particle  density  is  low  enough  (or  if  the 
electronic  charge  £  is  zero)  then  the  equation  reduces  to  the  usual 
simple  diffusion  equation. 


A  solution  of  Equation  3  was  desired  for  an  instantaneous  line 
source  in  an  infinite  medium  at  time  t  =  0.  It  can  be  shown  that  for  r 

very  small  a  solution  of  (3)  is: 

-r2/4  x  t 

_  N  e _ 

n  4  TOO  t  (4) 

where  N  =  the  number  of  ions/cm  introduced  by  the  line  source  and: 

N  e  2 

“  -  D(1+  —  >  (5) 

This  solution  approaches  the  usual  line  Source  diffusion  solution  if 
N  e  /kT  approaches  zero.  The  magnitude  of  N  £  /k  T  is  a  measure 
of  the  effect  of  electrical  forces  in  magnifying  the  diffusion  constant  D. 


An  approximate  solution,  applicable  to  all  values  of  r  and  t  is: 

2 


n  = 


N 


1 


m 


m  +  1 


(it  ,  12 _  •  •  •  -X-  \ 

l+x  +  ~  T  > 

4  7T  D  t  2  m. 


-x 


(6) 


x  = 


ft 


N  <£  A  T 


4  Dt 


The  differential  equation  does  not  impose  any  restriction  upon  the  choice 
of  m,  however  if  m  is  chosen  as 

m  =  N  e  2 A  T  (7) 

then  the  solution  at  r  =  0  is  the  same  as  that  given  by  (4)  and  (5). 

Further,  if  N  — ^-0,  the  usual  diffusion  solution  is  obtained.  In  general, 


B-5 


2 

N  €  /kT  will  not  be  an  integer  and  the  last  term  of  the  polynomial 
should  be  chosen  with  m  equal  to  that  integer  closest  to  Nfi*7"/  k  T. 


The  solution  expressed  by  (6)  is  very  flat  out  to  x  m. 

In  the  interior  region,  the  space  charge  effect  has  produced  a  nearly 
uniform  particle  density.  The  point  at  which  the  diffusion  front  is  the 
steepest  is: 


r 


2 


2(1  + 


2  Nfc 
kT 


Dt 


(8) 


Although  this  is  one  way  of  specifying  the  extent  of  the  diffusion,  an 
alternative  way  is  by  computing  the  mean  square  diffusion  distance: 


72“  =  2(2+  )  Dt  (9) 

Equations  (8)  and  (9)  are  to  be  regarded  as  useful  approximations.  They 

2 

yield  the  same  result  only  if  N  6  /  k  T  =1.  In  the  event  that 
2 

N  €  /  k  T  <  1,  the  spreading  is  controlled  primarily  by  diffusion  and  (9) 

is  a  natural  extension  of  the  usual  diffusion  distance,  r^  =  4  D  t.  If 
2 

N  €  /  k  T  >  1,  then  the  diffusion-like  solution  becomes  subordinate  to 

the  space  charge  effect  and  (8)  gives  a  better  physical  representation  of  the 
radius  which  contains  a  substantial  fraction  of  the  particles.  In  this 
case  the  distribution  is  fairly  uniform  over  most  of  the  confining  radflis. 
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NOTE  2: 


Mean  Distance  for  Electron  Attachment  to  Halogens. 


4/ 

The  data  of  R.H.  Healey  indicates  an  attachment  cross-section 
-18  2 

of  about  10  cm  for  1  volt. electrons  in  cl  2  anc*  ®r2  anc*  PerhaPs  3  times 

this  value  for  the  attachment  of  low  energy  electrons  (say  0.5  ev  or  less) 

to  I2.  The  data  does  not  clearly  indicate  what  cross-sections  might  obtain 

for  electron  energies  at  lower  energies.  It  will  be  assumed  that  the  attach- 

”18  2 

ment  cross-section  remains  at  10  cm  for  electrons  down  to  zero  energy, 

that  the  electrons  which  leave  the  filament  are  rapidly  scattered  by  an  inert 

gas  to  thermal  temperature  (300°K)  where  the  mean  electron  speed  is 
7 

1.08  x  10  cm/second. 

The  flux  of  electrons  is  <p  -  nv  and  <p  satisfies  the  diffusion 
equation."^ 


Here  \  is  the  scattering  mean  free  path,  is  the  absorption 

s 

(attachment)  microscopic  cross-section,  S  is  the  source  strength  and 
v  is  the  electron  speed.  Assuming  that  the  gas  is  largely  xenon,  then 
the  electron  scattering  cross-section  in  xenon  at  1  atmosphere  is  not 

less  than  3800  cmVctn^.^  Hence,  the  scattering  mean  free  path 

-4 

is  not  greater  than  1/3800  =  2.63  x  10  cm.  If  the  total  gas  is  made 
up  of  xenon  containing  a  small  fraction,  k,  of  halogen,  then  the 
absorption  cross-section  is: 

=  (2.44  x  1019)  (  10“18)  k 

=  24.4  k  cm'Vcm8. 

19  o 

(2.44  x  10 _ =  gas  particle  density  at  1  atmosphere  and  300  K). 

4/ 

As  reported  by  S.C.  Brown,  loc.  cit,  pg.  178 

5/  Glasstone  and  Edlund  "Elements  of  Nuclear  Reactor  Theory",  Van  Nostrand, 
19  52  ,  p.  101,  Equation  5.35.1. 

6  / 

See  S.C.  Brown,  loc.  cit,  Figure  1.10  and  the  remarks  following 
Equation  1.4. 
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We  are  primarily  concerned  with  the  mean  distance  from  the 
origin  at  which  electrons  become  attached  and  so  consider  the  steady 
state  situation  described  by: 


^r 


2 


+ 


d(f> 

d  r 


X — $  -  0 


This  equation  has  the  solution: 

<f)  =  &  KhQ  (  p  r) 

where  KhQ  is  the  Bessel  function  i  Hq  ^  (i  x)  tabulated  by  Janke  and 

Emde,  Fourth  Edition,  page  236.  The  function  is  positive  everywhere, 

infinite  at  the  origin,  and  decreases  in  exponential-like  fashion  with 

increasing  radius .  The  factor  p  =  V  3  /  \  =  528  ^^k  with 

the  previously  given  values  of  and  /\  .  The  solution  is  normalized 

s 

by  requiring  that  the  integral  of  21*  (fa  over  all  space  (i.e.  the  total  rate 
of  electron  attachment) ,  is  equal  to  the  rate  at  which  electrons  are  admitted 
at  the  source . 

N  =  d  v  =  ZTTCC'£aJ^  KhQ  (  ^  r)  r  dr 


Integrating  as  indicated, 

a  =  N/  4X>a=  (3/4)  N/  \  s 

The  normalized  solution  is  then: 

=  f  ■f’  Kho  (  ?  r) 

s 

Although  some  electrons  continue  as  r  becomes  large,  the  decay  of 

electron  flux  is  very  rapid.  The  number  of  electrons  absorbed  interior  to 

r  is:  ,r 

n  =  271  |  (p  r  dr 

' o 
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TC  N 
2 

Therefore:  n_  _ 

N 

For  n/N  =  0.9  (i.e.,  90%  of  the  electrons  have  been  absorbed),  ^  r  must 

satisfy 


[  £  -  V  Khi  <  *  r) 

1  -  Y  P  a,  (  ?  r) 


^  r  Khj  (  ^  r)  =  (0.1)  (2/n>  )  =  0.0636 

or  ^  r  =3.22 

with  (5  jt  528  V  k,  the  radius  r  within  which  90%  of  the  electrons 
become  attached  is  given  by  r  =  3.22/528  y~k.  This  radius  is  shown 

as  a  function  of  k  in  the  following  table: 


k 

0.01 

0.02 

0.04 

0.10 


r  (cm) 
0.0610 
0.0431 
0.0305 
0.0194 
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APPENDIX  C 


STUDIES  PERFORMED  ON  THE  THERMAL  SYSTEM 

1.  The  thermal  system  consists  of  a  pulsed  heater  and  a  down¬ 
stream  detector  which  senses  the  arrival  of  pulses  of  heated  gas.  The 
detector  may  be  a  thermocouple  or  a  resistance  thermometer,  the  heater 
is  a  simple  wire.  The  heater  and  sensors  must  be  small  enough  to  have 
little  thermal  inertia  at  the  operating  frequency  (30  cps).  The  transit 
time  of  the  gas  pulse  must  be  sufficiently  short  that  the  gas  pulse  does 

not  seriously  decay  by  thermal  conduction  between  creation  and  detection. 

2.  The  calculations  cf  this  appendix  show  how  the  important  dimensions 
(transit  time  and  wire  size)  are  established. 

3.  The  problem  of  the  behavior  of  the  heat  pulse  as  it  transits  from 
the  heater  to  the  detector  is  equivalent  to  the  problem  of  a  heat  source 

1/ 

moving  through  a  fixed  fluid.  This  problem  is  treated  in  the  literature 
and  the  equation  for  the  temperature  of  the  fluid  may  be  written: 

— “■  =  OcV 2  Q  S  (  x  -  vt)  •  £  (y)  (1) 

Specifically,  this  relates  to  a  line  source  of  power  Q  (t)  per  unit  length 
which  moves  at  speed  v  through  a  fixed  medium.  The  temperature  at 
x,  y,  and  t  is  ©,  oC  is  the  thermal  diffusivity,  p  c  is  the  volumetric 
heat  capacity  and  the  delta  functions  indicate  that  the  source  is  always 


^  See  Sneddon  "Fourier  Transforms"  McGraw-Hill,  1951,  p.  201 
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In  the  present  problem,  we  are  more  concerned  with  the  solution  as  one 
follows  the  source  at  a  fixed  distance  ^  =  v  t  -  x.  The  The  integral  is 
approximated  by  the  method  of  stationary  phase,  so  an  approximate,  algebraic 
solution  for  the  case  with  y  =  0  is: 

6  ^ >  -  b  pc  h* f-v -  •  [  2u)(  VJTi- »]  (3> 

and,  if  y  ^  0,  then 

at#  ‘  <*> 

The  nature  of  the  approximation  is  such  that  the  solution  is  only  valid  in 
the  region  behind  the  source  and  only  if  the  driving  frequency,  f  =  to  /2TC 
satisfies: 


In  the  present  case,  the  value  of  the  right-hand  side  of  (5)  is  about 
38  cps,  whereas  the  contemplated  operating  frequency  is  20  to  30  cps. 
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4.  The  solution  (3)  shows  that  at  distance  ^  from  the  source, 
the  temperature  oscillates  between  zero  and  peak  values  given  by: 


4  f  c  V  *  0(7  T 

'  2 
The  thermal  diffusivity 06  ,  is  taken  as  0.093  cm  /sec; the  volumetric 

-4  3  o 

heat  capacity,  c,  is  1.33  x  10  cal/cm  -  K,  and  the  speed  v  is 

16.35  cm/sec.  In  order  to  obtain  a  peak  temperature  rise  of  l°C  at  ^ 

0.5  cm  from  the  source,  the  peak  source  strength  (from  Equation  6  must 
-4  -3 

be  8.2  x  10  cal/sec.  or  3.42  x  10  watts  per  cm  of  heater  length. 
With  this  power  input,  the  peak  temperatures  observed  at  other  points 
^  ,  y  are  shown  in  the  table  below: 


© 


max 


)  = 


(6) 


9  (  *  )  °C 


f  =  0.32 _ (L40 _ 0.50 


0 

1.25 

1.12 

1.00 

0.02 

1.18 

1.07 

0.97 

0.04 

0.99 

0.93 

0.86 

0.06 

0.75 

0.74 

0.72 

(The  dimensions  |  and  y  are  centimeters) . 

5.  The  frequency  restriction  implied  by  (5)  is  somewhat  arbitrary 
since  it  does  not  quantitatively  specify  the  errors  arising  when  the 
equality  obtains.  For  this  reason,  numerical  integrations  were  performed 
using  an  operating  frequency  of  25.5  cps.  The  temperature  fluctuation 
(peak-to-trough)  at  ^  =  0.32  and  y  =  0  was  found  to  be  1.19  degrees 
instead  of  the  1.25  degrees  of  the  preceding  table.  This  indicates  that 
the  approximate  solution  is  quite  good  in  the  region  of  interest 
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(20  to  30  cps).  It  should  be  noted  that  if  the  heater  driving  voltage  is 
alternating  at  25  cps,  then  the  actual  heat  input  is  varying  at  a  frequency 
of  50  cps.  For  this  reason,  the  criterion  (5)  may  be  interpreted  as 
indicating  that  serious  attenuation  of  signal  frequencies  should  not  be 
expected  below  38  cps.  The  numerical  integration  indicates  that  signal 
frequencies  up  to  51  cps  are  propagated  without  serious  attenuation  for 
the  contemplated  system  parameters . 

Although  the  foregoing  establishes  the  fact  that  low  pressure 
gases  can  be  satisfactory  from  the  standpoint  of  thermal  conductivity,  it  has 
been  found  that  they  are  unsatisfactory  from  a  heat  transfer  view  point. 
Specifically,  they  do  not  remove  heat  rapidly  enough  from  the  heater 
win  and  so  lead  to  an  undesirably  long  thermal  pulse  (  <  50  milliseconds). 
Calculations  indicate  that  liquids  and/or  high  density  gases  should  be 
satisfactory  in  both  respects  and  current  activities  are  devoted  to  invest¬ 
igations  of  liquid  systems. 
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